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Voyager 1 and 2 data reveals that magnetic field fluctuations are compressive and exhibit a 



Gaussian distribution in the compressed heliosheath plasma, whereas they follow a lognormal 

distribution in a nearly incompressible supersonic solar wind plasma. To describe the evolution 

of magnetic field, we develop a nonlinear simulation model of a partially ionized plasma based on 

^ two dimensional time-dependent multifluid model. Our model self-consistently describes solar wind 

1*3 | 

plasma ions, electrons, neutrals and pickup ions. It is found from our simulations that the magnetic 
"c/3 ■ 
£>"V field evolution is governed by mode conversion process that leads to the suppression of vortical 

modes, whereas the compressive modes are amplified. An implication of the mode conversion 

process is to quench the Alfvenic interactions associated with the vortical motions. Consequently 

anisotropic cascades are reduced. This is accompanied by the amplification of compressional modes 

o' 

that tend to isotropize the plasma fluctuations and lead to a Gaussian distribution of the magnetic 
I/-} . field. 
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I. INTRODUCTION 

Space plasma is in a fully developed turbulent state |l|, |2j. Turbulent interactions are 
mediated by the solar wind that emanates from the Sun and propagates outwardly It 
interacts with partially ionized interstellar gas predominantly via charge exchange, and 
creates pick up ions |4j, |5j. Near the termination shock (which is about 90-100AU from the 
Sun), the supersonic solar wind decelerates, heats up, and it is compressed. It becomes 
subsonic in a region called heliosheath. In the heliosheath region, the solar wind plasma is 
compressed. The solar wind further interacts with interstellar neutrals via charge exchange. 
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These interactions are described comprehensively by Zank in Ref. |4J. During its journey 
from the Sun, the solar wind plasma develops multitude of length and time scales that 
interact with the partially ionized interstellar gas and nonlinear structures develop in a 
complex manner. Many features of the in situ heliosheath plasma have been surprising 
and were not expected from the existing analytic and simulation modeling. One of the 
most notable Voyager observations is the solar wind plasma near the heliosheath is subsonic 
and compressive |6J. The subsonic and compressed solar wind plasma exhibits a Gaussian 
distribution in magnetic field fluctuations contrary to the lognormal that is typically observed 
in the non compressive solar wind plasma [6| . The physical processes leading to the Gaussian 
distribution in magnetic field fluctuations are not understood. 

A primary goal of this paper is to describe a self-consistent evolution of the compressed 
solar wind plasma fluctuations by examining why magnetic field fluctuations exhibit a Gaus- 
sian distribution. For this purpose, we develop a fully self-consistent description of plasma- 
neutral coupled system and investigate compressive and non-compressive characteristic of 
magnetic field fluctuations in the context of partially ionized solar wind plasma. This issue 
is critically important in space plasmas because of its ramifications on origin of cosmic rays, 
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energetic particles, partially ionized turbulence and many other |4|, [5|, |7|. 

In section 2, we describe our new multi fluid model of plasma that is coupled with neutral 
gas in a partially ionized environment. Our model self-consistently describes the evolution 
of solar wind ions, electrons, pickup ions and neutral fluids. Implicit in our model is the 
interaction of small scale turbulence with a compressive plasma. Section 3 describes our 
simulation results dealing with the compressive characteristic of the solar wind plasma. Sec- 
tion 4 describes statistics of magnetic field fluctuations in compressive and non-compressive 



MHD plasma and finally section 5 summarizes our major findings. 

II. MULTIFLUID TURBULENCE MODEL 

Our nonlinear simulation model employs the dominant components of multi fluid species 
of the solar wind plasma. It includes plasma electrons, pickup ions, solar wind ions, and 
neutral gas. The solar wind ions interact with the interstellar neutral hydrogen via charge 
exchange that depends on the relative speeds of the solar wind and neutral atoms (^, l5|, LZj] . 
We assume that fluctuations in the plasma and neutral fluids are isotropic, homogeneous, 
thermally equilibrated and turbulent. The characteristic turbulent correlation length-scales 
(A c ~ l/k c ) are typically bigger than charge-exchange mean free path lengths (A ce ~ l/k ce ) 
in the space plasma flows, i.e A c ^> A ce or k ce /k c 3> 1. 

The fluid model describing nonlinear turbulent processes in the interstellar medium, in the 
presence of charge exchange, can be cast into plasma density (p p ), velocity (Up), magnetic 
field (B), pressure (P p ) components according to the conservative form 
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The above set of plasma equations is supplimented by V • B = and is coupled self- 
consistently to the neutral density (p„), velocity (V n ) and pressure (P n ) through a set of 
hydrodynamic fluid equations, 
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Equations ([TJ) to ([2D form an entirely self-consistent description of the coupled 
plasma-neutral turbulent fluid. The charge-exchange momentum sources in the plasma 
and the neutral fluids, i.e. Eqs. (TlJ and (|2]), are described respectively by terms 
Q M (U p ,V n ,pp,p n ,T„,Tp) and Q M (V n , U p , p p , p n ,T n ,T p ). A swapping of the plasma and 
the neutral fluid velocities in this representation corresponds, for instance, to momentum 
changes (i.e. gain or loss) in the plasma fluid as a result of charge exchange with the 
neutral atoms (i.e. Qm (U p , V n , p p , p n ,T n ,T p ) in Eq. (TlJ). Similarly, momentum change 
in the neutral fluid by virtue of charge exchange with the plasma ions is indicated by 
Qm(V„, XJ p ,p p , p n , T n , T p ) in Eq. (T5]). In the absence of charge exchange interactions, the 
plasma and the neutral fluid are de-coupled trivially and behave as ideal fluids. While the 
charge-exchange interactions modify the momentum and the energy of plasma and the neu- 
tral fluids, they conserve density in both the fluids (since we neglect photoionization and 
recombination). Nonetheless, the volume integrated energy and the density of the entire 
coupled system will remain conserved in a statistical manner. The conservation processes 
can however be altered dramatically in the presence of any external forces. These can in- 
clude large-scale random driving of turbulence due to any external forces or instabilities, 
supernova explosions, stellar winds, etc. Finally, the magnetic field evolution is governed 
by the usual induction equation, i.e. Eq. ([TJ, that obeys the frozen- in- field theorem unless 
some nonlinear dissipative mechanism introduces small-scale damping. 

Our model equations can be non-dimensionalized straightforwardly using a typical scale- 
length (£ Q ), density (p ) and velocity (t> ). The normalized plasma density, velocity, en- 
ergy and the magnetic field are respectively; p p = p p /po,Up = \J p /vq,P p = P p /pqVq^ = 
B/^oy/po. The corresponding neutral fluid quantities are p n = p n /po, U n = U n /t>o, Pn — 
Pu/poVq- The momentum and the energy charge-exchange terms, in the normalized form, 
are respectively Q m = Qm^o/Po^o; Q e = Qe^o/Po v o- The non-dimensional temporal and 
spatial length-scales are t = tv /£ , x = x/£ . Note that we have removed bars from the 



set of normalized coupled model equations ([I]) & (J2J). The charge-exchange cross-section 
parameter (a), which does not appear directly in the above set of equations, is normalized 
as a = n £ a, where the factor n £ has dimension of (area) -1 . By defining n ,£ through 
a ce = I/uqIq = k^ e , we see that there exists a charge exchange mode (k ce ) associated with 
the coupled plasma-neutral turbulent system. For a characteristic density, this corresponds 
physically to an area defined by the charge exchange mode 
expressions for charge exchange sources are taken from Refs 

III. SIMULATION RESULTS 
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We have developed a two-dimensional (2D) nonlinear fluid code to numerically integrate 
Eqs. flTJ to (|2j). The spatial discretization in our code uses a discrete Fourier representation 
of turbulent fluctuations based on a pseudospectral method, while we use a Runge Kutta 
4 method for the temporal integration. All the fluctuations are initialized isotropically 
with random phases and amplitudes in Fourier space. A mean ambient magnetic field is 
assumed to be present to describe the large scale background magnetic field in the plasma. 
This algorithm ensures conservation of total energy and mean fluid density per unit time 
in the absence of charge exchange and external random forcing. Additionally, V • B = is 
satisfied at each time step. Our code is massively parallelized using Message Passing Interface 
(MPI) libraries to facilitate higher resolution. The initial isotropic turbulent spectrum of 
fluctuations is chosen to be close to k~ 2 with random phases in all three directions. The 
choice of such (or even a flatter than -2) spectrum does not influence the dynamical evolution 
as the final state in our simulations progresses towards fully developed turbulence. While 
the turbulence code is evolved with time steps resolved self-consistently by the coupled fluid 
motions, the nonlinear interaction time scales associated with the plasma 1/k • U p (k) and 
the neutral 1/k ■ V n (k) fluids can obviously be disparate. Accordingly, turbulent transport 
of energy in the plasma and the neutral fluids takes place on distinctively separate time 
scales. 

We now analyze statistics of magnetic field fluctuations in both compressive and non- 
compressive MHD plasma to describe the statistics of magnetic field fluctuations. It should 
be noted that the initial fluctuations in our simulations comprise both the vortical (i.e. 
irrotational motion of fluid flow) and compressional (due to the longitudinal flow motion) 
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FIG. 1: Quantitive evolution of compressive (k com = \/J2k |k ' U p | 2 / J2k |U P | 2 ) and vortical 



(k vor = \/J2k |k x UplVSfc | Up | 2 ) components of plasma flow. Here k and Nk are respectively 
the wave vector and total number of modes. The summation is carried over the entire turbulent 
spectrum. Initially, the vortical component is large. As time progresses, compressive component 
dominates over the vortical counter part in the heliosheath. 

components. Our previous work show that the vortical component of fluid flow dominates 
over the compressive component in a supersonic solar wind plasma [8j. We use this result 
as a basis to develop a self-consistent description of compressive plasma fluctuations. In the 
latter, the vortical motion is sustained predominantly by shear Alfvenic modes that govern 
nonlinear cascade in the solar wind plasma. The compressive modes, on the other hand, are 
composed of fast/slow modes. The former survives, whereas the latter decays in the solar 
wind [8]. By contrast, the Voyager's observations indicate that the solar wind becomes more 
compressive in the heliosheath plasma [6|. Hence the compressive component of the flow 
is expected to dominate the vortical component in the heliosheath plasma. To understand 
this apparent discrepancy between the compressive and vortical modes in the solar wind and 
heliosheath plasmas, we follow the evolution of the two components in our simulations by 
initializing the velocity field with a higher magnitude of vortical component. Our simulation 
results are shown in Fig. (JT]) for 512 2 modes in a two dimensional box of length 2-k x 2tt. The 
other parameters in our simulations are; charge exchange k/k ce ~ 0.1 — 0.01, fixed time step 
dt = 10~ 3 , and collision parameter v ~ 0.1 — 0.001. The background constant magnetic field 
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FIG. 2: The magnetic field fluctuations in the non-compressive MHD turbulence exhibits a lognor- 
mal distribution in our simulations. The magnetic field fluctuations are dominated by the vortical 
modes. 



B = 0.5. Our simulations are fully nonlinear because the ratio of the mean and fluctuating 
magnetic fields 5B/B ~ 1. 

As the evolution proceeds, nonlinear interactions quench the vortical component and am- 
plifies the compressive counterpart. Consequently, the latter grows, while the former decays 
eventually and stays constant. The compression of the velocity field corresponds essentially 
to the compression of the magnetic field by virtue of the field and flow that are coupled 
strongly under the ideal frozen-in-field state. Our multi fluid simulations thus demonstrate 
that progressive development of compressive turbulence plays a catalyzing role in the mode 
coversion (vortical to compressive) process. Once the mode conversion process is over, 
the two components decouple permanently and evolve independent of each other. Further 
growth of the compressive component in our simulation is ascribed to turbulent fluctuations 
that are converted predominantly into the compressive mode by nonlinear processes whereas 
minimal or almost no flux of energy is transmitted into the vortical motion. Our simulations 
results, describing the predominance of the compressive modes over the vortical as shown in 



Fig. (JTJ) , are qualitatively consistent with the Voyager's observations 
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IV. STATISTICS OF MAGNETIC FIELD FLUCTUATIONS 

We next analyze the magnetic field fluctuations in our simulations. The results of our 
simulations, shown in Figs (J2J) & (J3J), describe probability distribution function (PDF) of the 
magnetic field fluctuations respectively in the non- compressive and compressive regimes of 
solar wind MHD turbulence. The PDF of the magnetic field in the non-compressive plasma 
is consistent with the lognormal distribution. This is shown in Fig. fl2]). By contrast, the 
magnetic field in the compressive plasma follows a Gaussian PDF as shown in Fig. (j3]). The 
latter is consistent with the Voyager 1 observations as reported by Burlaga et al. [6]. A 
Gaussian PDF corresponds typically to a uniform, random and isotropic distribution, and a 
mean deviation in any of the latter leads essentially to a skewed or lognormal distribution [9|. 
In the context of our simulations [see Figs ([1]) & (J2J)], we infer that a lognormal distribution 
of B in the non-compressive plasma results primarily by the predominance of vortical motion 
in magnetized plasma that primarily gives rise to Alfvenic-like fluctuations. In the presence 
of a mean or ba ckg round magnetic field, Alfvenic fluctuations tend to anisotropize the 



energy cascades 10|, llll]. Consequently, migration of turbulent energy is non symmetric 
along and across the mean magnetic field. The anisotropic cascade is therefore a process 
that could potentially lead to a skewed or lognormal distribution of the magnetic field in 
the non-compressive region which is dominated by nearly incompressible vortical motion. 
By contrast, plasma is dominated by the high frequency fluctuations in the compressed 
region. The effect of Alfven waves is relatively weak in this region as compared to the 
vortically dominated non-compressive plasma. Owing thus to the weaker Alfvenic effect, 
the compressional and relatively high frequency motions in plasma tend to isotropize the 
magnetic field fluctuations. Hence a Gaussian PDF follows in the compressive plasma. The 
physical process, describing how Gaussian and lognormal distributions occur respectively in 
compressive and non-compressive plasma, is consistent with our simulations shown in Figs 
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V. CONCLUSIONS 

In summary, we have investigated evolution of the magnetic field fluctuations in small 
scale compressive and non-compressive MHD plasma turbulence in a partially ionized en- 
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FIG. 3: Probability distribution function of the magnetic field in the post-shock is consistent 
with the Gaussian distribution. The magnetic field fluctuations are dominated by the compressive 
modes. 



vironment. Our results are useful in describing the magnetic field data from Voyagers 
We find that initial turbulent fluctuations, comprising both the vortical and compressive 
motion, evolve towards a state in which the vortical motion predominantly governs non- 
linear interaction in the non-compressive plasma by exciting Alfvenic modes. By contrast, 
the mode conversion process in the compressive plasma leads to the suppression of shear 
Alfvenic vortical modes whereas the compressive modes are amplified. The latter isotropizes 
the PDF of magnetic field in the compressive plasma. 

To summarize our findings, we find that the probability distribution function of mag- 
netic field in compressive MHD fluctuations is a Gaussian. The depleted vortical motions 
suppress the Alfvenic modes in the compressive MHD plasma. This we believe is one of 
the plaussible reasons why the magnetic field fluctuations are transformed into a Gaussian 
(from the lognormal) in partially ionized compressive solar wind plasma turbulence. Our 



fj| and theoretical predictions 
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results, consistent with the Voyager observations 

be useful in the context of heliospheric plasma where charge exchange interactions govern 

numerous features of the solar wind plasma J4J, |fj, [7]. 
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